Anarchy in the neutrino sector? 
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If the connection between the fundamental theory and neutrino mass matrices M v is sufficiently 
complicated all traces of underlying symmetries might be erased and an 'anarchical' M v , with 
random entries, can arise. It has been claimed that such random matrices generically prefer large 
mixing making unnecessary to invoke flavour symmetries to explain the observed pattern of neutrino 
mixing angles. A critical analysis of this idea (paying particular attention to the role of unphysical 
phases and coordinate dependence) shows that basis-independent, random M„'s are not biased 
towards large values of the neutrino mixing angles, contrary to previous claims. As a result, the 
case for anarchy in the neutrino sector is considerably weakened. 



PACS: 14.60.Pq, I2.60.Jv IFT-UAM/CSIC-03-18. 

1. Introduction. In a series of three interesting pa- 
pers [1-3] the idea that the measured neutrino parame- 
ters [4] could arise from a random mass matrix has been 
proposed and studied. A neutrino mass matrix with ran- 
dom entries could be the result of a fundamental theory 
that is itself complicated or that requires a convoluted 
mechanism {e.g. with several sources) for generating neu- 
trino masses. 

In ref. [1] , a random scan of 10 6 neutrino mass matrices 
(of Dirac, Majorana or see-saw type) is performed and 
the resulting distributions for mixing angles and for the 
ratio R = Amj? o; /Am^ tm are studied. The distributions 
for the atmospheric and solar angles peak at maximal 
mixing (which is welcome) while U e 3 is in general not 
particularly small (which can be a problem, or a promise 
for its measurement in the future). Finally, it is easy to 
get a small R (~ 0.1), especially in the see-saw case. 

Ref. [2] sheds lihgt on the previous results by realizing 
that the distributions of mixing angles are in fact dictated 
by the requirement of basis independence of the random 
scan. If M v is diagonalized by a mixing matrix U that 
belongs to a certain Lie group [e.g. U(S) for a complex 
3x3 Majorana M v ], then the probability distribution 
for the mixing angles that parametrize U follow the Haar 
measure of that Lie group. This nice property allows an 
analytic understanding of the predictions of anarchy. 

Finally, ref. [3] tests the hypothesis of anarchy in the 
neutrino sector by performing a Kolmogorov-Smirnov 
(KS) test taking as data the observed neutrino mixing 
angles (choosing the best fit points for the solar and at- 
mospheric mixing angles and leaving the CHOOZ angle 
free, but below its bound) . The test should measure how 
well the anarchy hypothesis describes the data or, alter- 
natively, see if the data are able to rule out this hypoth- 
esis (with some degree of C.L.). It is found that the KS 
probability in favor of the anarchy hypothesis is ~ 12%. 

In this letter, the predictions of the anarchy hypothesis 
for different types of neutrino mass matrices are carefully 
re-examined in section 2. This analysis shows that anar- 
chy has no preference for certain ranges of mixing angles. 
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In particular, there is no preference for maximal mixing. 
A critical analysis of statistical tests of the anarchy hy- 
pothesis, performed in section 3, confirms this result. 

2. Mixing angles predicted by Anarchy. One of 
the main ingredients of the hypothesis of anarchy in the 
neutrino sector is the requirement of basis-independence 
in flavour space: i.e. the entries of M„ are random num- 
bers distributed in some specified range according to a 
flat probability density, that remains flat in any basis. 
As an example, take a 2 x 2 real and symmetric M„. 
Factoring out the absolute neutrino mass scale m, take 
the entries of M u = M v /m as random numbers in the in- 
terval [— 1, 1] with uniform probability density. The dif- 
ferential probability density for each entry is simply eJMy 
and in the 3-D space {Mn, M12, M 22 } that parametrizes 
M v , the differential probability density is given by the 
3-form dP = dMu A dM 12 A dM 22 . It is straightfor- 
ward to prove the invariance of dP under rotations in 
flavour space. (This result can be extended to more 
general matrices: for complex entries one simply uses 
dMij = dSfte Mij A gKoti My.) A uniform probability 
density for the entries of M„ is not enough to guaran- 
tee basis-independence. The range in which these entries 
can vary must also be a basis-independent domain. In the 
previous example, this domain was a cube, which is not 
rotation-invariant, but one can simply take My inside a 
sphere of some specified radius: Tr MlM v < r 2 . 

Once the basis-independence of the random matrices 
is enforced, the mixing matrices U that diagonalize them 
are distributed in a well determined way: if the U's are 
matrices of a given Lie group, the differential probability 
density for U, d£2, is simply the invariant Haar measure of 
that Lie group (see [2] ) . In terms of the invariant Maurer- 
Cartan 1-forms uj a = -iTr T a XJ- x dXJ ", where T a are the 
generators of the Lie group, then dSl — e ai .a 2 ....a N ^ ai A 
... A uj aN , with TV the dimension of the group. 

Let us now examine in turn several examples of Majo- 
rana mass matrices (the results can be extended to Dirac 
or see-saw matrices) . 
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2.1 Real 2x2 M„ 

These matrices are diagonalized by rotations, U = 
R(0), so that the corresponding Lie group is U(l) and 
dfl = d6, which is explicitly rotation invariant. The dis- 
tribution of the mixing angle 9 is therefore flat and no 
mixing angle is preferred. The conclusions of ref. [2] are 
different: by writing d6 = dsin 2 2#/(2sin4#) it is shown 
that the distribution of sin 2 29 peaks at zero and maximal 
angles, which are claimed to be preferred. Such peaks, 
however, do not imply that angles near zero or maximal 
mixing are favored: these peaks are simply the result of 
using a particular coordinate in the circle, sin 2 29, which 
distorts the uniform distribution in angles. To see this 
even more clearly, notice that we have not specified in 
what basis we are measuring 9. We would find peaks for 
sin 2 29 in any basis, but if they were to imply a prefer- 
rence for maximal and zero mixing, one arrives at the 
paradoxical conclusion that different physical angles are 
preferred depending on the basis used. 

Given the flat distribution in 9, all one can say is that 
large mixing angles are not unlikely. The relative like- 
lihood of finding 9 in a region near zero mixing or in a 
region near maximal mixing is simply given by the rela- 
tive size of the two regions. In section 3 we discuss the 
statistical tests that might be applied to data in order to 
see how well they support the anarchy hypothesis. 

2.2 Complex 2x2 M„ 

These Majorana matrices are diagonalized by U(2) 
matrices that can be written as U = exp [1(77/2 + 
u!a3)]R(9)exp[iipas\, where 77 and to are unphysical 
phases, 9 is the rotation angle and <f> a CP violating 
phase. The distributions of these quantities in the an- 
archy scenario follow from the U(2) Haar measure 

dQ, = ds 2 e A dr\ A duo A d<p , (1) 

(with s x = sinx). Rewriting ds 2 , = s 2 gd9, ref. [2] con- 
cludes that, in this case, maximal mixing is preferred. 
This conclusion, which seems now inescapable, leaves still 
the paradox that the preferred angle seems to depend on 
the basis chosen to measure 9 (in spite of the fact that 
dCl is indeed basis independent). One might argue that 
the measured 9 is defined in the flavour eigenstate ba- 
sis and therefore such basis should be used. However, 
one should be allowed to determine the probability dis- 
tributions and the most probable angle in a rotated basis, 
provided one rotates back later to the flavour eigenstate 
basis. The paradox is that, doing so, the final result for 
the expected 9 depends on the intermediate rotated basis 
used. 

In order to understand why, according to (1), 9 = tt/4 
seems to be preferred while 9 = is disfavoured, it is 
useful to inspect the explicit expressions for the invariant 
Maurer-Cartan forms: 

to 1 = dr\ , lo 2 = s 2 8C24>du) - s 2(j> d9 , , . 

lo 3 = s 2 eS24>d.Lo + c 2< f>d9 , u 4 = c 2 gduo + d(p , 

with dQ — lo 1 A lo 2 A lo 3 A lo 4 . We see that, for 9 — > one 
has uj 2 A lo 3 — ► 0, while this is not the case for finite val- 



ues of 9 due to the presence of duo terms. In other words, 
the subspace of parameter space corresponding to 9 = 
has one dimension less than that for 9 = 9q ^ and 
the extra volume associated with the latter case opens 
up in the direction of the unphysical parameter lo. How- 
ever, all matrices related to each other by changes in an 
unphysical parameter should be considered as physically 
indistinguishable and therefore the preference for maxi- 
mal mixings does not stand on firm ground. 

Keeping the unphysical phases 77, lo helps dealing with 
basis independence but introduces a bias in the distribu- 
tions of physical parameters. How do we get rid of this 
problem? We simply fix 77 = lo = 0, which we are free to 
do for any M v . The crucial point is that this condition 
is rotation invariant, as can be readily checked. After 
fixing this 'gauge', 9 transforms under a rotation by an- 
gle a simply by 9 — > 9 + a while the phase </> does not 
change (notice that in the presence of non-zero unphysi- 
cal phases the transformation properties under rotations 
are more complicated). Therefore, after 'gauge-fixing', 
the only volume element in the space {9, (/)} that is rota- 
tion invariant must be of the form [5] 

dVL gf = T{<p)d9 A d4> . (3) 

Because invariance under rotations is not enough to 
fix F{4>) in (3), further assumptions besides anarchy are 
necessary to make any prediction about the distribution 
of the phase 4>. This situation is similar to that of the 
distribution of mass eigenvalues [2] which are unknown 
functions of mass differences. In analogy with that case, 
one might choose these functions in such a way as to get 
a volume element in parameter space that is simple when 
expressed in terms of mass matrix elements. Such choice 
is useful to generate samples of random matrices like in 
[1] but now in a 'fixed gauge'; however, its theoretical 
motivation is not clear. A simple choice for a rotation 
invariant volume clement (in terms of matrix entries) is 
dVgj = dKeAfiiAd9?eM 2 2AdKeMi2Adlog3mMi2, which 
leads to 

dVgf = -2(mi - m 2 )dmi A dm 2 A c 2 d</> A d9 . (4) 

This privileges small values of 4> but, as explained, other 
choices of invariant volume are possible. The uncertainty 
in the distribution of <f> does not affect the distribution of 
9, which is again flat: all mixing angles are equally prob- 
able and the comments at the end of section 2.1 apply 
also to this case. 
2.3 Real 3x3 M„ 

A symmetric real M v can be diagonalized by an orthog- 
onal matrix U belonging to 50(3). The parametrization 
U = R 2 3(9 1 )R 31 (9 2 )Ri2(93), with R i3 {9) a rotation in 
the plane ij by an angle 9, is quite convenient to discuss 
neutrino oscillations. The angle #3 is the relevant angle 
for the oscillations of solar neutrinos, 9 2 is the small mix- 
ing angle bounded by CHOOZ experiment and 9\ is the 
relevant angle for the oscillations of atmospheric neutri- 
nos [4]. The differential probability distributions for the 
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mixing angles 9i derived from anarchy is then dictated 
by the SO (3) invariant Haar measure 

dQ = d9 l A cos 6 2 d6 2 A d9 3 . (5) 

In order to understand better what are the implications 
of (5) concerning the likelihood of different values for the 
mixing angles, it is convenient to switch to a different 
parametrization of U. Let ft be the axis of rotation in 
flavour space [parametrized in terms of polar angles a, (3 
as ft = (spc a , sps a , cp)\, and 9 the rotation angle. We can 
then write U = R 12 (a)R 31 (f3)R 12 (9)R 31 (-f3)R 12 (-a). 
In terms of these parameters one obtains 

dfl = da A sm. f3dj3 Asm 2 (6/ 2)d6 . (6) 

For a fixed value of 9 this measure is just that on the 
sphere (da A sin (3d(3) with polar coordinates a, (3. Al- 
though this measure is zero for sin (3 = 0, this does not 
mean that rotations along the z axis are disfavoured but 
rather that the distribution of the rotation axis ft is uni- 
form in flavour space, as it should be the case if one starts 
with a rotation-invariant mass matrix. A useful analogy 
is to consider a uniform distribution of points over the 
surface of the Earth. In terms of longitude (a) and lat- 
itude ((3) the probability density has the same form as 
above and the expected number of points inside a cer- 
tain region is simply proportional to its area no matter 
where this region is on the surface of the Earth. 

The previous discussion shows that in order to make 
meaningful statements about what angles are preferred, 
one should take into account the area (or volume) of the 
regions of parameter space that are compared. What 
anarchy implies, as encoded in (5), is that all regions 
of parameter space arc equally likely The 'preference' 
for small 9 2 in (5) is due to the fact that the volume of 
parameter space with small 9 2 is larger [in the same sense 
that (3 ~ 0,7r/2 is disfavored in (6) because the area of 
the poles is smaller than the rest of the Earth's area] . 

2.4 Complex 3x3 M v 

A complex symmetric M v is diagonalized by a t/(3) 
matrix that may be written in this way: U = 
diag(e M Se M2 ,e^) R 23 (6i)e- iSX ' R 13 (6 2 ) e lSX ' R l2 (6 3 ) 
diagte-^+He* 6 ,^ 2 ) with A' = diag(l, 0, -l)/2. In 
this notation 011,2,3 are unphysical phases. The distribu- 
tion of angles and phases according to the original idea 
of anarchy would then follow the distribution 

dfl = dsl A dc 2 A ds\ Ad8 A d^i A dS, 2 A dai A da 2 A da 3 

which corresponds to the U(3) Haar measure. This re- 
sults from combining 9 invariant Maurer-Cartan forms 
(straightforward to compute but lengthy). 

As was discussed for the complex 2x2 case, one cannot 
use dfl above to extract physical implications without 
getting rid of the unphysical degrees of freedom (which 
otherwise affect the distribution of physical parameters). 
This problem is solved by 'choosing a gauge' like a\ — 
a 2 + 5 and a 2 — a 3 = — £ 2 , which is invariant under 



rotations of basis. The invariant volume element after 
this gauge-fixing is of the form [6] 

dQ gf = T(8 - 3(, 2 )g{s 5 s 3 c 3 ){c 2 d9 1 A d9 2 A d9 3 ) 

A(d& Ad£ 2 AdS/ss) ■ (7) 

This volume element is composed of three parts: T and 
Q are unknown functions of their invariant arguments; 
the angular element is the volume clement for real 3x3 
matrices (5), which corresponds to a homogeneous dis- 
tribution of mixing angles; and finally there is a phase 
volume element. The predictions of anarchy for the an- 
gles 9\ and 9 2 are similar to those of the real 3x3 case: 
all values of 9\ are equally likely; small 9 2 is preferred due 
to volume effects. Without further assumptions on the 
form of T and G, anarchy does not make more detailed 
predictions for 6*3 and the physical phases (except for the 
prediction of a flat distribution for £1). If Q(x) = x, dCl g f 
gives a preference for large solar mixing. 

3. Statistical test of the hypothesis of anarchy 
In order to assess how likely it is that the observed 
neutrino angles arise from a distribution dictated by the 
hypothesis of anarchy, ref. [3] performs an statistical test 
as follows. Consider for simplicity the solar angle 9 3 and 
assume that the probability density is flat in x = sin 2 9 3 
(as was the prediction for 3 x 3 complex Majorana ma- 
trices in ref. [3], revised in this paper) and take x as a 
random variable in the interval [0, 1] (the same analy- 
sis applies to a flat distribution in d9, with x — 29 /n). 
Experimentally, the best fit point for solar oscillations 
corresponds to x e = 0.3. Based on this, ref. [3] tests 
the anarchy hypothesis using standard methods, by com- 
paring the distribution function predicted by anarchy, 
F(x) — J Q f(x')dx' = x, with the 'empirical' distri- 
bution function, F e (x), obtained from the best guess 
probability density, f e (x') = 5(x' — x e ), that leads to 
F e (x) — Jq f e (x')dx' = Q(x — x e ). The mismatch be- 
tween the two distribution functions is measured by the 
(two-sided) Kolmogorov-Smirnov statistic 

D = max x [\F e (x) - F(x)\] . (8) 

The higher D is, the worst F(x) accomodates the ob- 
served data. More precisely, the probability of obtaining 
a worse value of D from a different observed x is 

where D e is the KS statistic for x e . It is clear that, 
according to this procedure, the value x e = 1/2 (i.e. 
maximal mixing) is the best possible result in favor of 
the anarchy hypothesis. 

This is a puzzling result. How can a single data point 
be able to provide evidence in favor of a flat distribution 
at all? Why should data values close to x e = 1/2 be 
the best support for a flat probability distribution for xl 
While this seems reasonable for x — sin 2 9 3 it is certainly 
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not for x = 29/w, and ref. [3] uses the above test, with- 
out modification, also for distributions flat in x = cos 9. 
There are three problems with the previous analysis. 

The first is that x is an angular variable (0 < 9 < n/2 is 
mapped into the interval [0, 1] with x = and x = 1 cor- 
responding to the same physical point) and the test used 
is suited for linear variables. In fact the statistical analy- 
sis of angular or directional data requires some departure 
from the usual linear statistical methods and there is a 
well developed corpus of knowledge on this topic (for an 
introduction to this subject see e.g. [7]) which is relevant 
in many areas of science, from the analysis of directional 
patterns in bird migration to magnetic pole wandering, 
coincidence of planetary orbital planes or distribution of 
earthquakes in the globe, among others. The version of 
the Kolmogorov statistic that is best suited to angular 
random variables is Kuiper's test [7], defined as [8] 

D = max [F e {9) - F{9)] + max e [F{9) - F e {9)] . (10) 

This statistic has the desirable property of being inde- 
pendent of the choice of origin for 9 (see Appendix for a 
simple proof), which is not the case for the linear Kol- 
mogorov statistic, Eq. (8). 

However, the previous improvement is not sufficient in 
view of the second problem: the smallness of the sample. 
One can easily see that the Kuiper statistic for a single 
random draw, 9 e , is D e — 1, which is independent of x e . 
Although this seems to be a reasonable result (it implies 
that a single data point cannot give evidence in favor or 
against the hypothesis of a flat probability distribution) 
it cannot be taken seriously because the Kuiper test is 
always D e = 1 for a single data point irrespective of the 
distribution being tested! One possible way out would 
be to improve the test by taking into account the differ- 
ence between the distribution functions over the whole 
interval ('integral' variants of the Kuiper test exist [9]) 
but, as we show below, these improvements would also 
be insufficient. 

The third problem is that the K-S test is invariant 
under reparametrizations of x, e.g. it cannot tell apart 
distributions that are flat in x = sin 2 9 from those that 
are flat in x = 29 /ir, and this is inappropriate. 

An alternative approach (for x = 29 /n) is to enlarge 
the data sample by considering the origin, 9 a = 0, de- 
fined by the flavour basis, as another data point. In that 
case the probability of getting a fit worse than the ob- 
served one, according to the Kuiper test, is given again 
by Eq. (9) with x e = 2(9 e — 9 )/tt, which is explicitly 
rotation invariant. From this result one would conclude 
that 9 e — 6*o ~ 7r/4 would give the strongest support to a 
flat distribution and 9 e — 9 ~ 0,7r/2 the weakest, which 
agrees with the conclusions of Ref. [3] . However, this re- 
sult cannot be trusted. The reason is that even two points 
are too few to give evidence in favor of a flat distribution 
function: any set of two data points has the same proba- 
bility of any other (this situation changes only for three 
or more data points). 



We therefore abandon the attempt at performing a 
meaningful statistical test of the anarchy hypothesis, at 
least for uniformly distributed angles (in principle a par- 
tial test could be devised for angles which are not uni- 
formly distributed, like #2) and prefer to limit ourselves 
to the following estimate, similar to the analysis per- 
formed in Ref. [1]. Assuming a large number of random 
neutrino mass matrices, we can estimate now what frac- 
tion of them will have angles in some neighbourhood of 
the observed values (say in a subspace v of the total vol- 
ume V). For 3x3 real Majorana matrices [or for complex 
ones with Q{x) = 1 in Eq. (7)], that fraction will be given 
simply by the ratio J v dQ/ Jydil. Using ranges similar to 
those in Ref. [1] (sin 2 29 2 < 0.15, sin 2 26> 1:3 > 0.5) the in- 
teresting area is about a 4.8% of the total. This number 
is similar to the one obtained in Ref. [1] for the real case 
(in which case we both agree on the angle distributions) 
before applying the cut in mass splittings. If Q(x) = x 
is chosen in (7), then the preference for maximal solar 
mixing raises the previous percentage to 6.8%. Still, we 
regard this as a very small basis upon which to build the 
case in favour of anarchy. In contrast, attempts based on 
flavour symmetries fare considerably better [10,11] and 
are physically much more appealing: they imply that 
what we learn at low-energy from measurements of the 
parameters in the neutrino sector will be a window (with 
a view) to fundamental physics at much higher energies. 

Appendix To see that the Kuiper's test is invariant 
under rotations, note first that, due to the 'periodicity' 
property of angular distribution functions, [F(9 + 2tt) — 
F(9) + 1], any angular interval [a, a + 2ir] is suitable 
to compute D in Eq.(10). If we now change from the 
random variable 9 to the rotated 9' = 9 + a, the new 
distribution functions are F'(9') = F(9 + a) - F(a) and 

D 1 = max e ,[F' e {6') - + max e ,[F'(9') - F' e (9')] 

= max [F e (9 + a) - F(9 + a)} - F e (a) + F{a) 
+ max e [F(9 + a) - F e (9 + a)} + F e (a) - F(a) 
= max e [F e (9) - F(9)} + max e [F(9) - F e (9)} , 

so that D' = D. A less straightforward proof of this 
property can be found in [7] . 
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